Free energy barrier for single-chain melting and crystallization 
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In this paper, we report dynamic Monte Carlo simulations of melting and crystallization in a 
single-chain system. Their free energy barriers are calculated by the umbrella sampling method and 
can be described well by a simple expression AF = nAf + a(N — n) 2 ^ 3 , where n is the amount 
of molten bonds, Af is the free energy change of each molten bond from a crystalline state, TV 
is the chain length, and a is the surface free energy of crystallite. We found that, together with 
the expression AF = nAf + a(N — n) 1 / 2 for molecular nucleation, the molecular-weight dependent 
properties of the free-energy barriers for polymer primary and secondary nucleation, in particular, 
the molecular segregation during crystal growth, can be reproduced. Then for the mechanism of 
polymer crystallization, we suggested a quantitative model of intramolecular nucleation, as a direct 
development from the previous qualitative description of molecular nucleation model. 



I. INTRODUCTION 

Polymer nucleation dominates the crystallization pro- 
cess in a quiescent liquid. jl| The nucleation rate G can 
be described as the formula proposed by Becker: Q 

G = G exp[(AE + AF c )/(k B T)}, (1) 

where Go is a prefactor, AE is the activation energy for 
diffusion across the phase boundary, AF C is the free en- 
ergy barrier for critical nucleus formation, and ks Boltz- 
mann constant, T temperature. After the primary nu- 
cleation, the secondary nucleation is believed to be a 
rate-determined step of polymer crystal growth in the 
bulk states, leading to a constant linear crystal growth 
rate and its proportionality with exp(— 1/T/AT). ||] 
For small molecules, the secondary nucleation is a two- 
dimensional nucleation on a smooth growth front as a 
key process to generate layer by layer advance. The well- 
known Lauritzen-Hoffmann theory for polymer crystal 
growth was suggested based upon this kind of secondary 
nucleation process. Q However, the large-scale smooth 
front may not be available due to the thermodynamic 
roughening at high temperatures, and the ill conforma- 
tion of chains on the growth front may construct an en- 
tropic barrier for polymer crystal growth rather than the 
surface free energy penalty. Moreover, the qualitative 
model of molecular nucleation was suggested to explain 
the molecular segregation with different chain length dur- 
ing polymer crystal growth. Q 

Nevertheless, most of molecular-weight-dependent 
properties of polymer crystallization have not been ex- 
plained well by these incompatible models. Long chains 
require for large supercoolings in both primary nucle- 
ation |^| and crystal growth than short chains, while the 
superheating in melting seems quite small and insensi- 
tive to the chain length. The linear crystal growth 
rate of polymers has been found to be significantly chain- 
length dependent. [00] However recently, a master curve 



for the temperature dependence of crystal growth rate 
was found, which implies that the chain-length depen- 
dence should be attributed to the prefactor Go in Eqn. 
[j] rather than the free-energy barrier AF C for secondary 
nucleation. |l(J This chain-length independence was also 
observed in the primary nucleation rates of polyethy- 
lene, (ll]] The complicated molecular-weight-dependent 
properties imply a different crystallization mechanism of 
polymers from that of small molecules. Therefore, under- 
standing these properties is of fundamental importance. 

In principle, the molecular-weight-dependent proper- 
ties can be attributed to the processes of diffusion or the 
free energy change, which are sensitive to the size of the 
whole molecules. These kinds of process can be the in- 
tramolecular process since all the units of a molecule can- 
not crystallize simultaneously, especially for long flexible 
chains. When a polymer system performing melting and 
crystallization is limited to a single chain, the free en- 
ergy change of intramolecular process becomes essential. 
A crystallite formed solely by a single chain has been 
observed in experimen ts l|12| and reproduced in many 
numerical simulations. |13[| Recently, benefited from the 
simple lattice model [ ff4[ and the umbrella sampling tech- 
niques Jig] , we can calculate the free energy change for 
phase transitions of a single-chain system through dy- 
namic Monte Carlo simulations. []l6|| In this paper, the 
chain-length dependence of the free energy barrier for 
phase transitions was investigated. 

The rest of this paper is organized as follows. Af- 
ter the simulation technique section, a simple estima- 
tion of the free energy change for a single-chain sys- 
tem was performed. Then, the results of estimation for 
the barrier of single-chain equilibrium states and of pri- 
mary nucleation were compared to the simulation re- 
sults and experiments, and the formula was extended 
to a secondary nucleation process. Most of important 
molecular-weight-dependent properties of polymer crys- 
tallization and melting can thus be discussed. After that, 
we suggested an intramolecular nucleation model for the 
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mechanism of polymer crystallization and discussed its 
advantages. 



II. SIMULATION TECHNIQUES 

In our simulations, the lattice chain occupying a num- 
ber of consecutive sites performs micro-relaxation in a 
cubic box with periodic boundary conditions. |fL7f The 
size of box is 256 x 256 x 256, which is large enough to 
avoid the contact from one side of the suspending single 
coil to the other side through the periodic boundaries. 
The self- and mutual-avoiding jumps of the monomers on 
the lattice can be along either grid lines or diagonals, so 
the coordination number of cubic lattice is 9 + 8 + 9 = 26. 
Only the parallel packing tendency E p for each pair of the 
non-bonded bonds in neighbor is employed as the driv- 
ing force for crystallization, while the bending penalty 
along the chain and the mixing heat of the chain units in 
a solvent are set to zero for simplicity. [|l4| The system 
temperature is then represented as a dimensionless term 
T* = ksT/Ep in the Metropolis sampling. 

A series of single-chain systems containing variable 
chain units (16, 32, 64, 128, 256, 512, 1024) were in- 
vestigated. The sample, prepared initially in a regularly 
folded state, relaxed down to the equilibrium state at the 
preset temperature. The degree of melting is traced as 
the number of molten units, which is defined as the bond 
containing less then five bonds packing around in paral- 
lel. The free energy distribution to the number of molten 
units was estimated according to the Boltzmann relation 



AF/{k B T) = -lnP(n), 



(2) 



where P{n) is the probability to find a sample having 
n molten units. The umbrella sampling method permits 
us to calculate the free energy distribution in a window 
away from the most favorable state in equilibrium. 1 18| A 
parabolic biasing potential was employed here. [jl9| The 
results for free energy distribution in several such kind 
of windows covering over a wide range of n were then 
merged together through parallel tempering. |2C|] 



III. THEORETICAL ESTIMATION 

The free energy change of a single-chain system can be 
estimated with a fully ordered ground state in the bulk. 
In the simple lattice model, the potential energy loss can 
be estimated as AE = nE p (q — 2)/2, where n is the num- 
ber of molten units, q is the coordination number of a reg- 
ular lattice, the first 2 is the number of connective bonds 
along the chain, so q — 2 is the possible number of the 
parallel packing bonds around each bond, and the sec- 
ond 2 is a symmetrical factor for bond-bond pair interac- 
tions. The interfacial free energy penalty is estimated to 



be <t(N — n) 2 / 3 , where N is the total number of bonds on 
the chain, and a is the surface free energy. The entropy 
gain should be attributed to the conformational entropy 
change of the molten bonds, and as a first approximation, 
their conformation can be in analogy with a self- avoiding 
walk. [^IJ Therefore, AS — nlnu + (7 — l)lnn. Here, \i is 
the connective constant, which is the number of possible 
directions for each step of random walk, and for chain 
conformation, [i = q — 1; 7 is a critical exponent. The 
term Inn is a correction term for the self-avoiding walks 
deviated from the random walk and be actually equiva- 
lent to the mixing entropy of a single coil. This nonlinear 
Inn term is expected to be much smaller than that linear 
n term when n is large. Although the molten units may 
belong to several segments of the chain tied on the sur- 
face of crystallite, their linear terms contributing to the 
conformational entropy are additive. Therefore, we ne- 
glected the nonlinear term for simplicity. The free energy 
change of a single-chain system can be estimated as 



AF = ^-^nE p + a(N - n) 2/3 - k B Tn\n(q - 1) 

= nA/ + (r(A-n) 2 / 3 , (3) 



where Af = E p (q- 2) /2 - k B T\n(q - 1) is the body free 
energy change of each molten unit. Corresponding to our 
simulation model, N is fixed, q — 26, only a and T are 
required to be determined. 

The formula of Eqn. || reflects the Gibbs phenomeno- 
logical expression for nucleation of phase transitions. At 
the critical size (top of free energy curve), the free en- 
ergy change can be calculated from dAF/dn = as 
NAf + 4cr 3 /(27A/ e 2 g ). At the fully crystalline state, 
n = and hence the free energy change is erA 2 / 3 , while 
at the coil state, n = N and hence the free energy change 
becomes NAf. Therefore, the height of free energy bar- 
rier can be calculated from the difference of free energy 
changes between the initial and the critical states. 

Before starting the next section, several terms should 
be clarified as schematically shown in Fig. [l| When the 
lowest point of free energy change at the left-hand side 
(the crystalline state) has the same level as that at the 
right-hand side (the melt state), the sample system is 
in equilibrium between melting and crystallization. This 
temperature is the equilibrium melting point T eq . The 
free energy barrier at this temperature is the equilibrium 
free energy barrier AF eq . At a slightly high temperature 
T, the crystalline state becomes metastable, and the free 
energy barrier for melting is AF m with a superheating 
AT m — T—T eq ; at a slightly low temperature T, the melt 
state becomes metastable, and the free energy barrier for 
crystallization is AF C with a supercooling AT C — T eq — T. 
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FIG. 1. Schematic showing free energy curves in the vicin- 
ity of equilibrium melting temperature T eq , the critical free 
energy barrier AF eq , and the free energy barriers for melting 
AF m with a superheating AT m and for crystallization AF C 
with a supercooling AT C . 



IV. RESULTS AND DISCUSSION 

1. Single-chain primary nucleation 

Figure 2 shows the free energy changes for phase tran- 
sitions of a single 1024-mer, starting from the initial crys- 
talline state located at zero-line level. The equilibrium 
melting point was found to be 2.967Ep/ks- Deviating 
from this temperature, there will be a supercooling or a 
superheating to drive crystallization or melting respec- 
tively. 
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FIG. 2. Free energy curves vs. the number of molten units 
for 1024-mer at the denoted temperatures with the unit of 
Ep/ks- The numerical results (solid lines) are calculated by 
umbrella sampling with 15 windows. The dashed line is calcu- 
lated from Eqn. H with a = 15 E v and the fitted equilibrium 
melting temperature 3.2657-Ep/fcs. 

The height of equilibrium free energy barrier in simula- 
tions can be fitted by the simple estimation of free energy 
change in Eqn. ||. As demonstrated by the dashed curve 



in Fig. g, a proper set of a and T can fit the height 
of free energy barrier well. The significant deviation at 
the small amount of molten units should be attributed 
to some molten units of simulations located at the fold- 
end surface of the crystallite and the nonlinear entropy 
contribution in theoretical estimation neglected for sim- 
plicity. We are not going to improve the precision of 
estimation at small molten units, since here we only fo- 
cus our attention on the height of free energy barrier for 
phase transitions. 

At the equilibrium melting temperature, the coex- 
istence of the crystalline and amorphous states gives 
criV 2 / 3 = NAf eg , or Af eq = crA^ 1 / 3 . The height of 
equilibrium free energy barrier can be estimated as 

^ eq ~ 27Ap eq ~ 27 • ( > 

So the height of equilibrium free energy barrier is ex- 
pected to increase monotonically with the chain length. 
This is exactly what we observed in Fig. where the 
simulation results for the height of equilibrium free en- 
ergy barrier with variable chain length are quite close 
to the fitting curve of Eqn. || with just a single set of 
parameter a = 15E P . 
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FIG. 3. Height of equilibrium free energy barrier vs. chain 
length. The solid curve is calculated from Eqn. ^ with 
a = 15E P . 

The monotonic chain-length dependence of equilibrium 
free energy barrier implies that if the chain length is small 
enough, the height of this barrier can be smaller than the 
thermal fluctuation level. In other words, short chains 
should show no hysteresis in the temperature scanning 
for melting and crystallization. This can be identified 
in our simulations as demonstrated in Fig. |], where in 
short-chain cases, the phase transitions are reversible and 
appears continuous with temperature changes. Actually, 
this continuous change should be attributed to the aver- 
age of packing potential energy in each step of temper- 
ature change. If we trace an isothermal process in the 
vicinity of phase transition temperature of 64-mer, as 
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shown in Fig. [^, the single-chain system is quickly jump- 
ing between two extreme states and shows all-or-none 
feature. So during phase transition of short chains, the 
temperature change just shifts the relative probability of 
two states from one to another and leads to a continuous 
change in the average of packing potential energy. 



AF m = A/A - crN 
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FIG. 4. Heating and cooling curves (dashed and solid lines 
respectively) of the potential energy with the denoted vari- 
able chain length (number of units). The potential energy 
is denned as the mean value of 24 — p, where p is the num- 
ber of parallel bonds packing around each bond. The heating 
and cooling programs are the steps in the value of _E p /(fcsT), 
each having a length of 0.01 and a period of 10 6 Monte Carlo 
(MC) cycles. One MC cycle is defined as one jump for each 
chain-unit on average. The reported data are averaged over 
1000 samples, each with 1000 MC-cycle interval. 
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FIG. 5. Isothermal curve of the potential energy for 64-mer 
at T = 1.845i5p/fc_B. The initial state was picked up from the 
corresponding heating curve in Fig | after 10 6 MC-cycle an- 
nealing at this temperature step. 



When the temperature leaves away from the equilib- 
rium point, the height of free energy barrier for melting 
can be estimated as 



4cr 3 
27A/2 



(5) 



meanwhile, the height of free energy barrier for crystal- 
lization is 



AFr = 



4(7 J 



(6) 



27A/ 2 ■ 

At a fixed temperature (and hence A/), AF m depends 
upon the chain length, but AF C seems not. 

The chain-length independence of the free energy bar- 
rier for primary nucleation is in consistence with re- 
cent experimental reports about polyethylene primary 
nucleation rates. [jll| This reveals the importance of in- 
tramolecular process in polymer primary nucleation. Ac- 
tually, this independence is also verified in Fig. |[ where 
the cooling curves show the crystallization temperatures, 
which are insensitive to the chain length when the chain 
length is large enough to show the hysteresis. In the free 
energy calculation of simulations, we fixed the crystal- 
lization temperature there and found that the height of 
free energy barriers for crystallization is really insensitive 
to the chain length, as demonstrated in Fig. [| 
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FIG. 6. The free energy changes of single-chain systems 
with the denoted variable chain length (number of units) at 
the same temperature T = 2.174E p /kB- Curves are shifted 
to meet at their tops. 

The supercooling or the superheating required for 
phase transitions are the temperatures at which the ther- 
mal fluctuation can get over the free energy barrier for 
phase transitions and generate crystallization or melting 
spontaneously. The chain-length independence of AF C 
at the crystallization temperatures does not mean the 
same independence at the fixed supercoolings, since the 
equilibrium melting point, i.e. the reference temperature 
for supercoolings, may be sensitive to the chain length. 
Figure ^a shows the estimated free energy barriers for 
crystallization vs supercoolings, as well as that for melt- 
ing vs superheatings according to Eqn. [5] and |(| At a 
supposed thermal fluctuation level shown in the dotted 
curves, the height of free energy barrier for crystallization 
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is much more sensitive to the chain length than that for 
melting. This is also consistent with the experiments for 
the bulk polymer systems, which show that long chains 
require large supercoolings for primary nucleation, while 
the superheating required for melting seems insensitive to 
the chain length. Therefore, the rising of supercool- 
ing for bulk primary nucleation of long n-alkanes can be 
associated to their intramolecular chain-folding process. 
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FIG. 7. Height of free energy barriers for primary nucle- 
ation or melting (a) and secondary nucleation or melting (b) 
(AF C vs. AT C (solid lines) and AF m vs. AT m (dashed lines) 
with denoted chain length (number of units), calculated from 
Eqn. | and |(a), and Eq. | and | (b) with a = 15E P . The 
dotted line is the supposed thermal fluctuation level. 

2. Extension to secondary nucleation 

The success of the simple free energy estimation for 
polymer primary nucleation encourages us to extend our 
free energy estimation to secondary nucleation. Follow- 
ing the idea of molecular nucleation, the first part of 
molecule growing onto a smooth front will incur a penalty 
for a two-dimensional nucleation. So corresponding to 
Eqn. |, the free energy change of a single-chain growth 
onto a smooth front is 

AF = nAf + a(N -n) 1/2 . (7) 

And correspondingly, the height of free energy barrier for 



melting can be estimated as 

AF m = NAf - trN 1 ' 2 



4A/' 



(8) 



and the height of free energy barrier for crystal growth is 

„2 



AF r = 



4A/ 



(9) 



Again, at a fixed temperature (and hence A/), AF m de- 
pends upon the chain length, but AF C seems not. The 
chain- length independence of free-energy barrier for crys- 
tal growth agrees well with the experiments as mentioned 
master curves in the introduction section, ficfl Similar to 
primary nucleation and melting, the chain-length depen- 
dence of supercooling and superheating required to gen- 
erate crystal growth and melting, see Fig. [?]b, is also in 
consistence with experiments. 
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FIG. 8. Height of free energy barriers for crystal growth as 
a master function of temperatures estimated from Eqn. ^| and 
Eqn. [l^. The dashed line is indicating a fixed crystallization 
temperature for discussion. 

The equilibrium barrier AF e q can also be obtained 
with the same formula as AF r , 



Aft, 



4A/, 



(10) 



eq 



which implies a master curve for their temperature de- 
pendence. As demonstrated in Fig. |§|, if the sample sys- 
tem is fixed at a temperature indicated by a dashed ver- 
tical line, corresponding to an equilibrium melting point 
at this temperature, there is a critical chain-length N c for 
a polymer mixture with wide chain-length distributions. 
For the polymers with chain length larger than N c , they 
have the same height of free energy barrier for crystal 
growth at this temperature (see Eqn. ^), while their bar- 
rier for melting is much higher. This means that they can 
grow onto the front with enough stability. On the other 
hand, for those polymers with chain length smaller than 
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N c , they cannot stay on the growth front because of their 
lower barrier for melting than that for crystallization. So, 
N c is a critical value for molecular fractionation during 
polymer crystal growth. According to the equilibrium be- 
tween two states at this temperature, Af eq N = aN 1 / 2 , 
we obtain that T eq N^ 1 / 2 . This linear relationship can 
be well verified by the plot of experimental data for the 
equilibrium critical molecular weight of polyethylene melt 
fractionation at a crystallization temperature, see Fig. 0. 




FIG. 9. Equilibrium critical molecular weig ht W~ 1/2 of 
polyethylene melt as a function of crystallization tempera- 
ture T c . The data points are from Ref. 3 page 101 in Table 
V-7, Type M series, excluding the low molecular-weight end. 
The fitting line shows a correlation coefficient of 0.9997. 

Figure || also demonstrates the molecular segregation 
in the crystallites grown at different stages of a cool- 
ing process. When the sample system is cooled from 
high temperatures, the critical chain length is shifting 
to smaller. In a polymer mixture with a wide chain- 
length distribution, the longest chains will first approach 
the critical chain length on cooling, and therefore, con- 
tribute to the thick dominant lamellae, whereas those 
short chains approach the critical chain-length at the 
late stage of cooling and then contribute to the thin sub- 
sidiary lamellae. 

3. Intramolecular nucleation model 
Our estimation about free energy change of single- 
chain melting and crystallization gives a perfect quan- 
titative description to the molecular nucleation model. 
However, the original molecular nucleation model sug- 
gested only for the first part of molecule that incurs a free 
energy barrier for crystal growth onto the front. Actually 
in polymer crystallization, after the first part of molecule 
has finished the nucleation process, the fast growth of 
the rest part may quite easily be stopped, such as by the 
entanglement with other chains or by the impingement 
of the two-dimensional nucleus with other nuclei on the 
front surface. The entanglement has been found to make 
long-term pauses during single-chain crystal growth in 
the semi-dilute solution. |22] The impingement can be 



associated to the granular texture on the growth front 
observed recently under atomic force microscopy (AFM). 

n 

Therefore, the rest part of chain may still need to 
perform secondary nucleation process for further crys- 
tal growth. If this rest part of chain is longer than the 
critical chain length, as schematically shown in Fig. [ll], 
the free energy barrier for secondary nucleation, which is 
independent of the rest chain length, can still be got over 
in subsequent crystal growth. Since the coil size of the 
long-chain molecule can be much larger than the thick- 
ness of the lamellar crystallite, the multiple nucleation 
process may not be necessary to happen one-by-one in 
sequence like shown in Fig. 10, rather, they can occur 



simultaneously when the nuclei are in enough distance on 
the same front or belong to different crystallite fronts. 
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FIG. 10. Schematic of the multiple nucleation processes of 
single-chain crystal growth. The dashed vertical lines are indi- 
cating the stops of crystal growth due to some kinds of reason 
such as chain entanglement in the melt or nuclei impingement 
on the front. 

Based upon the above considerations, we suggest an in- 
tramolecular nucleation model. According to this model, 
any part of a molecule, if it is longer than the critical 
chain-length, can perform secondary nucleation process. 
If the rest part of chain is shorter than the critical chain- 
length, it will become a cilium tied on the fold surface or 
lateral surface of the lamellar crystallites. This multiple 
nucleation model can reproduce a reasonable trajectory 
of a single chain in a semi-crystalline contexture, where 
one molecule may pass through several lamellae. 



V. CONCLUSION 

Simple free energy estimation for the barriers of single- 
chain melting and crystallization can fit well to our simu- 
lation calculations and be consistent with the experimen- 
tal observations about their chain-length dependence. 
Then, most of important molecular-weight dependent 
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properties about polymer melting and crystallization be- 
come rational under the consideration of intramolecu- 
lar nucleation model. This proves the importance of 
intramolecular process in both primary nucleation and 
secondary nucleation processes of polymers. The advan- 
tages of intramolecular nucleation model, which directly 
develops from the idea of molecular nucleation model, 
include 1. reproduce the temperature dependence of nu- 
cleation rate as InG 1 /T/AT 2 for primary nucleation and 
InG 1/T/AT for secondary nucleation; 2. in a compar- 
ison with Lauritzen-Hoffman model, avoid the assump- 
tions such as the barrier of the first stem and a large 
smooth front for nucleation of new layer-growth; 3. in a 
comparison with Sadler-Gilmer model, the detailed pro- 
cess of crystal growth is much clear; 4. reproduce molec- 
ular segregation during polymer crystal growth; 5. re- 
produce the chain-length dependence phenomena of su- 
percooling required for primary nucleation and crystal 
growth or superheating required for melting; 6. repro- 
duce the chain-length independence of free energy barrier 
for both primary nucleation and secondary nucleation; 7. 
the multiple nucleation process can generate a reasonable 
trajectory of a molecule in the semi-crystalline contex- 
ture; 8. the impingement of the intramolecular nuclei on 
the front may lead to a granular texture; 9. reveal the 
importance of chain entanglement on the crystal growth 
rate. 

The diffusion process may contribute to the molecular- 
weight dependence of nucleation and crystal growth rate. 
It includes both the chain motion in the crystalline re- 
gion and in the melt region that beyond the scope of our 
current discussion for free energy changes. The lamellar 
thickness is a result of late stages of crystal growth, so 
we did not discuss here. The intramolecular nucleation 
model does not reject the intermolecular nucleation pro- 
cess, which may be dominant for the crystallizaiton from 
the dilute solution, especially for those short chains. We 
expect to get much more insights about the mechanism 
of polymer crystallizaiton from the intramolecular nucle- 
ation model in near future. 
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